Fractional statistics of topological defects in graphene and related structures 
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We show that fractional charges bound to topological defects in the recently proposed time-reversal-invariant 
models on honeycomb and square lattices obey fractional statistics. The effective low-energy description is 
given in terms of a 'doubled' level-2 Chern-Simons field theory, which is parity and time-reversal invariant 
and implies two species of semions (particles with statistical angle ±7r/2) labeled by a new emergent quantum 
number that we identify as the fermion axial charge. 



Introduction - When the excitations of a many-body sys- 
tem carry electric charge that is smaller than the charge of 
its constituent particles (e.g., electrons) the charge is said to 
be fractionalized. This phenomenon is known to occur in 
the fractional quantum Hall (FQH) liquids quintessen- 
tial strongly correlated systems with broken time-reversal 
symmetry, where fractionally-charged excitations also obey 
fractional exchange statistics |2]. Recently, two model sys- 
tems have been introduced on the honeycomb and square 
lattices (|3l |4|] that exhibit charge fractionalization without 
breaking of the time-reversal symmetry. These models, in 
essence, generalize the concept of fractionalization in poly- 
acetylene f§\ to two dimensions and, remarkably, can be con- 
sidered weakly correlated. The experience with FQH sys- 
tems suggests that the exchange statistics of these fractionally 
charged excitations could be anomalous. This question is in- 
teresting for several reasons. First, a very general argument 
can be made |j^ that would seem to prohibit the existence 
of anyons in systems that obey time-reversal symmetry. Sec- 
ond, fractional statistics have recently captured attention due 
to their relevance to topologically protected quantum informa- 
tion processing |7]. Since the honeycomb lattice is found in 
natural graphene |8], and the square lattice model could be re- 
alized in artificially engineered structures |9], the possibility 
of realizing anyons in time-reversal invariant systems has both 
theoretical and practical significance. 

In this Letter we construct the low-energy effective theory 
for the fractional particles in models lHlil]. We find that they 
are indeed anyons, albeit of a very special kind, described by a 
doubled\J(l)2 x U(l)2 Chern-Simons (CS) theory previously 
discussed by Freedman et al. ifioll . In its topological sector the 
theory contains two species of semions, which transform into 
each other under parity and time reversal, thus escaping the 
constraints imposed by the argument of Ref.lll Systems under 
consideration here 13|,|4i] represent the first explicit example of 
models for which such a gauge structure emerges as the low- 
energy effective theory. 

Fractional charge - The low-energy theory for fermions on 
the graphene honeycomb lattice 1 1 1.1 and the square lattice 
threaded with jt flux per plaquette 11211 is the Dirac Lagrangian 



£^ {f/[i^+me-'^^')if/ 



(1) 



where i. - jfiZf^, /u = 0, 1, 2, are 4 x 4 Dirac matrices 
in the Weyl representation and is a four-component Dirac 
spinor whose components index two Dirac points ('valleys') 



and two sublattices. The mass m arises from the dimerization 
of hopping amplitudes introduced in Refs. '3', I4I The phase 
X describes the direction of the dimerization pattern. Hou et 
al. jsll made a remarkable discovery that a topological defect 
(a vortex in;^') binds fractional charge ±e/2. Determination 
of this fractional charge relies on establishing the existence of 
unpaired zero modes in the solution of the associated Dirac 
equation For our discussion it is useful to de- 

duce the fractional charge by a method that does not rely on 
zero modes but instead exploits the long-distance behavior of 
Dirac fermions in topologically non-trivial backgrounds. 

The idea here, originally due to Goldstone and 
Wilczek fl6^, is to follow the flow of charge as we in- 
troduce a vortex into x by adiabatically deforming the 
order parameter, starting from a uniform configuration, 
e.g., = 0. This adiabatic insertion is made possible 

by temporarily enlarging the parameter space of masses by 
adding to the Lagrangian ([T]i a "y3-mass" term, m^iff'/itf- 
Physically this term corresponds to staggered on-site potential 
+m^ for the fermions on two sublattices of the square or 
honeycomb lattice. Such a term could in principle appear in 
the physical Lagrangian but we add it here by hand to enlarge 
the symmetry of the order parameter from U(l) to 0(3). Let 
us denote this 0(3) order parameter by the vector 



ip(x) — (m cos;^', m sinx, mj) 



(2) 



with a fixed length = - m 



- m^ and in the direction of 



the unit vector (fiix). Now we can adiabatically create a vortex 
in 01 +i(pi by first rotating (p from the initial uniform state "up" 
or "down" to (0, 0, sgn(m3)) and then flattening to the m^, - 
plane away from the center as illustrated in Fig. [Ua)- This 
creates a meron (half a skyrmion) in (p with the core pointing 
in the (0, 0, sgn(m3)) direction. On the lattice we can always 
take the limit of zero core size, thus recovering an ordinary 
U( 1 ) vortex. In the continuum description under consideration 
here, however, there always remains a single point (coincident 
with the singularity in x) where m^, retains a nonzero value, 
+M. This value distinguishes between the two diff'erent ways 
of creating a vortex and will be seen below to have physical 
implications. 

The utility of this procedure lies in the fact that the forma- 
tion of a vortex is completely smooth and we can thus calcu- 
late the charge accumulated during this adiabatic process by 
a perturbative loop-expansion of the fermion current opera- 
tor On symmetry and dimensional grounds we may expect 
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FIG. 1: (Color online) (a) Adiabatic formation of a vortex from a 
uniform, vortex-free configuration. The arrows represent the direc- 
tion of if) along the y = Q cut through the vortex center, (b) Feynman 
diagram used to calculate the fermion current. Solid lines represent 
the fermion propagator G{p) = - m)"'. 

that the fermion current will take the form of the conserved 
topological current in the 0(3) nonlinear cr model, 

yjT'' = -^e^vAeabclfiadv^bdAlfic, (3) 

up to an overall prefactor As in Ref.il6j the correctness of this 
anticipation and the prefactor can be established by comput- 
ing the expectation value of the fermion current = i/'Tf/'A in 
response to a small perturbation (p{x) imposed on top of a uni- 
form background (po, which we take to lie in the (1, 0, 0) direc- 
tion. A straightforward evaluation of the diagram in Fig.lTJb) 
yields (j^) = {M /^7:\M\^)ef^vAdv(fi2dAf3- Considering other 
permutations of external fields in the above diagram we iden- 
tify O/j) - f^'' ^'^^ general background field. 

To calculate the charge associated with the meron configu- 
ration let us substitute (p(x) from Eq. (|2]) into ([3]) to obtain 

0'/i> = ].., ^piy^(dym3)(dAX)- (4) 
47r|M| 

The total charge Q - J d^'^ijo) - 
(4-n\M\y^ J d^xm2(x)eijdidjX where we have used inte- 
gration by parts and the fact that m^ix) = far from the 
meron center. Recalling that for a static meron centered at the 
origin eijdidjxix) = 2jm5{x) with n the integer vorticity of x 
we obtain Q - insgn(m3(0)). This result has the satisfying 
property of not depending on the exact spatial profile of niT, 
but only on its asymptotic value at the vortex center. This is 
in complete agreement with Refs. lH, [l3l [3, when we 
recall that the charge of a n = +1 vortex can be either H-| or 
-J depending on whether or not is the zero mode occupied 
by an electron. In our present construction the sign of m-i is 
seen to play the role of the electron occupancy, despite the 
fact that the zero mode never enters our discussion. 

For future reference we also record the form of {j^) in the 
limit of infinitely small meron core, i.e., m^ix) - every- 
where except m3(0) = +M. Eq. (|4]i becomes 

iiii) = ^sgn(m3)e^v^5v5^. (5) 

This generalizes naturally to the case of multiple vortices with 
m3 - +M at each vortex center. 



Effective theory - We now turn to the discussion of statis- 
tics. Our strategy will be to find the low-energy effective the- 
ory for vortices in the Lagrangian (HJ. For the subsequent 
considerations it will be useful to modify the Lagrangian by 
including a coupling of the fermions to two gauge fields, 

X = lA ('^ - 4 + 75^ + me-'^^' ) i/r. (6) 

Afj and Bjj couple minimally to the electric and the axial 
fermion currents, = iAT/j^A ^nd y'g^ - i/'y^ygi/', respectively, 
and can be therefore identified as the ordinary electromagnetic 
field and the chiral gauge field introduced in Ref. 17. In what 
follows we treat A and B as static external fields that help us 
keep track of the charge content of various fields. 

We note that the coupling between vortices and fermions 
can be written as m{ij)+e'^4r- + (^_e"*i/'+) where if)± = i(l + 
75)1/' are the chiral components of the Dirac fermion. This 
suggests that in order to focus on the vortex degrees of free- 
dom we perform a singular gauge transformation t//-t — > 
e'^'^-t//±, where we have defined an (as yet) arbitrary partition- 
ing - X++X- of the vortex phase, such that upon encircling a 
vortex the spinor field remains single-valued lUal . (This sim- 
ply means that the phase field associated with any given vortex 
is assigned to either oi'X--) Under this transformation the 
Lagrangian becomes 

J:^^{W-i + yiH)+m)il/, (7) 

with shifted gauge fields 

a^i = A^ + ]jdf,{x+-X-)^ (8a) 

W = B^,+ ^-d^{x^+X-). (8b) 

A transformation of this type has been previously employed 
in the studies of the vortex state of a c/-wave superconduc- 
tor [fl8^. Its chief advantage is that it recasts a somewhat un- 
wieldy Lagrangian O) with a twisted mass term in the form of 
a gauge theory (|7]) which can be analyzed by standard field- 
theory methods. Specifically, our goal is to integrate out the 
fermi fields in (|7]l to arrive at the effective action in terms of 
the gauge fields which now encode the vortex degrees of free- 
dom. As is usually the case this procedure encounters ultra- 
violet divergences. These must be regularized in a manner 
that is consistent with the symmetries of the underlying lattice 
model. In the following analysis we shall pay particular at- 
tention to the time reversal and parity which we expect to be 
preserved in the low-energy theory. 

Symmetries - Besides Lorenz and gauge invariance, the 
low-energy theory also respects discrete symmetries. At m - 

parity, (f, x, y) (t, -x, y), acts as i/r ^ Pi// with P e 7iu(2) 
where u(2) is the group generated by {1, iyi,j5, 7375 ). Charge 

conjugation, C, can be similarly worked out to be — > Ct//* 
with C e 72u(2). The antiunitary time-reversal operation, 

T T 

(t,x,y) — > (-t,x,y), is given by — > Ti// with T € 7iu(2). 

s 

The Lagrangian is odd under the unitary operation (f, x, y) — > 
(-t,x,y) that sends iff ^ S iff with S E 7()U(2). That is, S 
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anticommutes with the Hamiltonian. This is identified on the 
lattice as a sublattice symmetry [3J which renders the energy 
spectrum symmetric around zero. 

r 

The sense of the vortex is switched under parity: x — > -X- 

9 

This means that {bQ,bi,b2) — » (-bo,bi,-b2)- The gauge 

field a must, however, behave as a vector like the usual elec- 

p 

tromagnetic field, (00,01,02) — > (oo, -01,02). So, we re- 

p c 
quire that x± ~X^- Similarly we find x± -X± ™d 

T T 

X± — > +X^- The latter implies (bQ,bi,b2) — > (-^0,^1,^2) and 

T 

(00,01,02) (00,-01,-02). These conditions also uniquely 
determine the operation of !P, C, T and S on the fermi fields 
in the presence of vortices (m + 0) to be 



P = nrs; c = 7275; ^ = nrs; s = 7073- 



(9) 



It follows that the bilinear 1^731^ is even under any operation, 
while \jjy^yf,\^ is even under S and odd under V, C and T. (In 
fact, the latter is true in the absence of vortices as well.) 

Topologicol terms - Armed with the above analysis we can 
now ask what topological terms are allowed by symmetry in 
the low-energy effective action for o and b. It is easy to see 
that conventional CS terms o-(dxo) and b-(dxb) break both 
P and T and are thus prohibited. One can, however, construct 
a mixed CS term. 



-£cs = -r'^ -(dxb) 
An 



(10) 



which obeys !P and T and is thus allowed. The value of the 
coeflicient k is tied to the value of fractional charge. The sim- 
plest way to see this is to note that in view of Eq. dSal i varying 
the action Sqs, - J dfd^xXcs with respect to gives the 
fermion current. 



6S, 



CS 



6A, 



^ S-^fivAdvdAX- (11) 
A,B=0 on 



Consistency with Eq. (|5]l then requires k = 2sgn(m3). 

The mixed CS term (fTOl l can also be obtained more directly 
from the Lagrangian ^ if we regularize it by adding a Pauli- 
Villars mass term. Based on our discussion of symmetries 
we choose to add the 73-mass term, which preserves both T 
and !P but breaks S. The physical amplitudes are found in a 
standard perturbative expansion using the Pauli-Villars sub- 
traction, :?lphys = ^Kg(m-i = 0) - y[reg('«3 — » °°). Integrating 
out the Dirac fermions to one-loop order is a lengthy but ulti- 
mately straightforward exercise that yields ifTgll 

r /a n2 \^\u2 Sgn(m3) 
Xeif = -7^(5 X of + —b^ + — O ■ (d X b). 

ll\m\ In In 

The last term is the mixed CS term obtained before. The 
(d X o)^ term describes the expected dielectric response of 
the Dirac medium. In the absence of the chiral gauge field 
the b^ term becomes simply (dxf. This reflects the cost of 
spatial and temporal variations in x(x) and encodes the usual 
logarithmic interaction between vortices. If the chiral gauge 
field is present and described at the bare level by a Maxwell 



term, then the interaction between vortices is exponentially 
screened at long distances, as in ordinary type-II supercon- 
ductors. Note that m appears as the charge of o and at the 
same time as the mass of b. Therefore, by tuning m — > we 
could expect to find a phase where o is massive and b is soft, 
hence the role of axial and regular currents is reversed. This 
would be a superconducting phase. 

Exchonge stotistics - We now analyze the implications of 
the mixed CS term ( fTOl i for the vortex statistics. To this end we 
substitute o and b from Eqs. dS) into Xcs and set A = B = 0, 



Xcs = (m+ ■ 5 X m+ — m_ ■ 5 X M_), 

\6n 



(12) 



where we defined - d^x±- It is now easy to understand 
the exchange statistics. Take two vortices at Xi and X2 and let 
the second one go on a path C2 around the first one, which 
remains static. Let us assume they both belong to the same 
±-partition and sgn(m3) > 0. In Eq. ( fT2] i we may then write 
u±-ui+ U2 and u^-^ where ^(5 x mj;)^ = (dx^/df)<5*^'(x - 
Xi:(f)) is the current density of the ^-th vortex. The topological 
phase e^'® that is accumulated in this process can be calculated 
from the u\ ■ dxu2 cross term in Eq. (fT2b : 



26* = +- f d^x5(2'(x-X2(f))— M 
2 J df 



= ±2 J* dx2 ■ ui = +n. 



I/' 



(13) 



This means that two such vortices behave as semions with the 
statistical angle - ±nj2. Alternatively, if we have two vor- 
tices in two different +-partitions there is no cross term and 
they will be mutual bosons. 

How can the assignment of vortices into the seemingly ar- 
bitrary +-partitions produce observable effects? The answer 
lies in the realization that this assignment actually entails a 
genuine physical distinction between vortices in the two par- 
titions. This can be seen by computing the axial current by 
varying >Scs with respect to B^, as in Eq. ( fTTI ). to obtain 
(js/i) = ^^jivxdvdx(x+ ~ X-)- This impHes that elementary 
vortices assigned to different partitions carry opposite axial 
charge Qs - ^ cfx{jso) - +k/4: If we had a chiral gauge 
field probe at our disposal, we could in principle detect the 
axial charge associated with a vortex just as we can detect its 
electric charge. In the absence of such a probe the axial charge 
still can affect the physics, e.g., by influencing the exchange 
statistics. 

The Lagrangian ([T]) posesses an emergent global symmetry 
under the transformation t// e'^^^^ij/ and;if X ~ ^Xo- This 
guarantees conservation of the axial charge in all low-energy 
processes. When a vortex is created, e.g., in a pair-creation 
process, it is endowed by a particular value of the electric 
charge (Q = +j) and the axial charge (^5 = +i). These val- 
ues then uniquely characterize the vortex; in particular they 
determine its exchange statistics. In order to better appreciate 
this key point it is useful to recall that ij/j- represent the field 
operators associated with the two different Dirac nodes [3, 4]. 
Physical meaning of the axial charge then becomes clear from 
the expressions jo - + lA-'A- and j^o - 'A+i/'+ - 
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Since the nodes interchange under T and V the axial charge 
is odd under these operations. The spatial components of 
are closely related to the 'valley currents' recently studied in 
the context of graphene 120*1. 

The above result (fTSl l could also be understood using the 
following physical picture. Starting with two vortices in the 
+-partition we could choose to transfer one, say vortex 2, to 
the — partition. The symmetric field, h, does not change but 
the anti-symmetric field, a, shifts by -ui. This shift could be 
absorbed in A, which will then attach a 2n flux to vortex 2. 
As a result, we pick up an extra Aharonov-Bohm phase tt by 
taking the flux 2ti of vortex 2 around the charge ^ of vortex 1 . 
The statistics remains unchanged. 

Doubled CS theory - It is possible to rewrite the Lagrangian 
(fT2l i in a more familiar form by introducing two auxiliary 
gauge fields Jl± mediating the statistical interaction between 
vortices, 

Xcs = ^l-^cry[^-5xy[^ + ;^-J?lJ (14) 

a-=± ^ ' 

where j± = j^(d x dx±) are currents associated with vortices 
in two partitions. Gaussian integration over leads back 
to Eq. (fT2l i. We recognize the above Lagrangian as the level-2 
doubled CS theory of Ref. 10, constructed there based on very 
general considerations. 

Concluding remarks - The statistical angle njl of a vacant 
vortex means that fusing two vacant vortices must result in 
a boson (with statistical angle 2n). At first sight this seems 
counterintuitive, since the resulting charge is | + j = +1 which 
could be expected to be a regular hole, and thus a fermion. 
But there is a subtlety here. Fusing two vortices gives a dou- 



ble vortex which supports two exact zero modes ^ 0, [Till . 
Now, suppose we remove the hole by adding an electron, so 
we have a neutral vortex. Since the electron can go to either 
of the two modes there is a pseudospin-^ degree of freedom 
attached to the neutral vortex. By the spin-statistics theorem 
this indeed should be a. fermion, consistent with the notion that 
the original double vortex is a charge- 1 boson. 

It is also interesting to contemplate exactly how our model 
evades the argument of Ref. 6. The crucial point is that since 
the axial charge is odd under time reversal the many-body 
ground state of the system with two vortices must be at least 
two-fold degenerate, with ^F* = and mutually orthogo- 
nal. However, the absence of anyons in a 7"-invariant system 
only follows from Ref. ^ when the state is non-degenerate. 
It is straightforward to generalize the argument to the dou- 
bly degenerate situation fl9^. One finds that because of the 
additional structure introduced by the degeneracy, 6 - +7r/2 
semions are allowed in addition to fermions and bosons, in 
agreement with the results of the effective theory. 

While this work was in the final stages a preprint by Cha- 
mon et al. ll2lll appeared, in which the 0(3) topological cur- 
rent is derived consistent with our Eq. Their conclusions 
regarding the exchange statistics, however, appear to disagree 
with ours. 
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